Estimate for the 0++ glueball mass in QCD 
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Abstract 

We obtain accurate result for the hghtest glueball mass of QCD in 3 
dimensions from lattice Hamiltonian field theory. Using the dimensional 
reduction argument, a good approximation for confining theories, we suggest 
that the 0+"^ glueball mass in 3+1 dimensional QCD be about 1.71 GeV. 
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The spectroscopy of QCD in the pure gauge sector, i.e., the glueball 
masses attracts considerable attention. In the quenched approximation, these 
glueball are non qq gluonic bound states formed by strong self-interactions 
of the gluons, and their masses vary from about 1.4 Gev to 2.5 Gev. The 
flux tube models, bag models, sum rules and lattice techniques have been 
used to extract the glueball masses, among which lattice QCD seems to 
give more reliable estimates. Experimentally, a lot of glueball candidates 
such as 6(1440), /o(1520), ^//j(1720) and ^(2230), produced in the J/i; 
radiative decays ||l|, |, H are within this range. The difficulty in experimental 
identification of a glueball comes from the complexity in determining the 
quantum numbers J^^ of these particles. 

The Monte Carlo simulation of lattice QCD on the Euclidean lattice has 
become a powerful and conventional method in spectrum calculations. Con- 
cerning the lightest glueball 0^"'', numerical data of 15 years ago, on rather 
small lattices with poor signal to noise ratio suggested M(0"^'^) ~ 1 Gev. The 
value of M(0"'"^) seems increasing with the lattice volume. Most recently, 
more accurate calculations by the IBM group on much larger lattices, 
higher statistics and better algorithm gave M(0^~'") ~ 1.740 ± 0.071 Gev, 
where the infinite volume extrapolation has been made. For a review of the 
current status, see ref. 0. 

Here we would discuss an alternative way 0, |^ to extract the glueball 
masses and wavefunctions by solving the Schrodinger equation p| 

H\F)=eF\F). (1) 

Whereas our goal is to do concrete computations 0] in four dimensional 
QCD, in this paper we would like to discuss quantitatively the properties of 
three dimensional QCD, an interesting and relevant but much simple theory. 



We will then use the idea of dimensional reduction [|1^, |TT|, to estimate 
the 0"''+ glueball mass in four dimensional QCD. 

With the standard notations {g being the lattice gauge coupling, a the 
lattice spacing, Ui = expiigaAi) the gauge link, Ei the color-electric field on 
the link /, Up the plaquette variable), the lattice Hamiltonian is given by 

^ = E ETET - ^ E Tr{Up + Ul - 2) (2) 

is the lattice version of the Yang-Mills Hamiltonian / d^^^x{e'^E ■ E + e^^B ■ 
B)/2 with e the continuum gauge coupling, E and B the color electric and 
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magnetic fields respectively. The glueball wavefunction in (|1|) is 

(3) 

created by the gluonic operator F{U) with the given quantum number J^"-^ 
acting on the vacuum The vacuum wavefunction \Q) satisfies H\Q) = 
e^\fl). An estimate of the glueball mass is then Mjpc = Ae = ep — e^. 

In a series of papers f^, 0, §], we developed a method for solving the 
lattice Schrodinger equation (|l|) in a new scheme which preserves the correct 
continuum behavior at any order of approximation. Physically, if one wants 
to well describe a glueball on the lattice, the size or the Compton length of 
a glueball, which is usually of the same order as that of a hadron, should be 
greater than the lattice spacing a. In other words, the low energy spectrum 
originates mainly from the long wavelength excitations. Our starting point 
is to embody this physical implication in the eigenvalue equation when solve 
it approximately. The philosophy is to have the correct long wavelength 
limit at any order of approximation. The advantage and reliability of such a 
method has been confirmed by the results of two-dimensional a models (in 
three-dimensional U(l) (in 0), SU(2) (in [§, |T|, 0]) and SU(3) (in 



[0) H) gauge theories: the results converge very rapidly, and even at very low 
truncation orders clear scaling windows for the vacuum wavefunction and 
mass gaps have been established. It is very exciting that for the a models, 
2-1-1 D U(l), and 2-1-1 D SU(2) gauge theories, the vacuum wavefunction 
and the mass gap are in perfect agreement with the most recent Monte Carlo 
data. In our pioneering study of 2-1-1 D SU(3), we obtained the first estimates 
for the vacuum wavefunction and the glueball masses P). Most of these 
results have been summarized in |jl^, p!8| . 

We begin with recapitulating briefly our method. Suppose the ground 
state has the form \Q) = exp[R{U)]\0) , with |0) being the fluxless bare vac- 
uum and R{U) being a linear combination of gauge invariant gluonic opera- 
tors G. Substituting this, (||) and (0) into we have the following exact 
eigenvalue equation for the operator F{U) 

J:{[Ei, [El, F{U)]] + 2[E,, F{U)] [E^, R{U)]} = ^-^F{U). (4) 
I 9 



In 0, H, |18|, we have illustrated how to obtain R{U) and F{U) by 
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expanding them in order of graphs (Wilson loops) Gn,i{U): 

F{U) = FniU) = fn.GnAU), (5) 
n n,i 

with n being the order of the graphs. In practice, equation (^) has to be 
truncated to some finite order N 

N 

I n ni+n2<N 



^^E^f/). (6) 

9 n 

from which the coefficients fn,i are determined. 

The first term in doesn't create higher order graphs, while the second 
term generates new or higher order graphs of order rii + n2. Therefore, one 
should carefully truncate the second term in this calculation. The essential 
feature of our approach is in the correct treatment of this second term. It has 
been generally proven p that in the long wavelength limit this term should 
behave as 

[Ei,F,{U)][Ei,R,{U)] oc a"^ Tr(P^^,,)2. (7) 

Not to violate this behavior, when the equation (^) is truncated to the Nth 
order, all the graphs created by [Ei, Fi{U)][Ei, Rj{U)] for + ^2 > must 
be discarded. For example, equation (j^) truncated to = 2 is 

[EuF, + F,]] + 2[EuF,][EuRi]} = ^(F^ + F,). (8) 
T ^ 

For N=3, the truncated equation (^) is 

Y{[Ei, [El, F, + F2 + Fs]] + 2[Ei, F,][Ei, R,] + 2[Ei, F,][Ei, R2] 
I 



+ 2[Ei, F2] [El, R,]} = ^(Fi + F2 + F3). (9) 
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Higher order truncated eigenvalue equations satisfy the same rule. From 
an eigenvalue equation at order N, we derive a set of nonlinear algebraic 
equations for the coefficients fn,i of the operators. Solving these equations, 
we obtain not only the wavefunction of the glueball, but also the glueball 
mass at order N. (For the vacuum wavefunction, see ref. 0). It is worth 
mentioning another advantage of such an approach: no group integration 
is necessary so that higher order calculations are feasible. The difference 
between different truncation orders is the estimate for the systematic error 
in the calculation. 

For a non-abelian gauge theory, the element A of the gauge group has 
to satisfy the uni-modular condition (or Caley Hamilton relation ||19||), 
which is for SU(3) 

Multiplying it by A^^, and then summing over the i index, it becomes 

Multiplying it again by epqr and summing it over the p index, we obtain 
AiiAkj = AijAki - TrA^{6ji6ik - SkjSu) - + - Aj^u + 4/'^u-(12) 

These formulae are useful in classification of graphics. Because of these 
conditions, one should choose properly an independent set from the graphs 
generated by the second term in (||). Mathematically, any independent set 
chosen in this way can be used in the calculation. 

Physically, the connected graphs represent more coherence and have less 



mixing with lower order graphs. It was shown in that the use of con 



nected set makes the convergence of the results much faster than the use of 



disconnected set in a (2+l)-dimensional SU(2) model. 

For the realistic gauge group SU(3), the complication is that not all the 
disconnected graphs can be transformed to the connected ones. However, 
we observed that if more disconnected graphs are transformed according to 
the uni-modular conditions ([T0|), (|lT]) or ([1^) into the connected ones [0], 
the scaling behavior was much better. We have also tested several sets of 
operators pO|. One of them is classified according to inverse of the graphs 



oc Inv[Gn,j] generated by the second term in (^ with the uni-modular 
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conditions taken into account. (For the definition of tlie inverse operator, see 



ref. 10 



Since QCD3 is a super-renormalizable gauge theory, the renormalization 
requirements amounts to dimensional analysis. In the weak coupling region 
(for large j3 = Q/g'^), because the renormalized charge e and the bare coupling 
are related by = e^a, dimensional analysis tells us that the dimensionless 
masses aMjpc should scale as 

aMjpc Mjpc _ , 

> — — ^ const., (13) 

g-^ 

from which the continuum physical glueball masses Mjpc are extracted. 

Using the techniques in p, ^ and after a careful analysis of our results, 
we obtain the value for M(0^"'")/e^ more accurate than the previous paper 
and estimate the systematic error due to the finte truncation to be 
less than 0.06. The validity of equation (|^) extends from /3 = 5 to /3 = 12, 
and in this range 

M(0++) , , 

\ ' ^ 2.15 ±0.06, (14) 

where the error denotes the systematic uncertainties due to the finite order 
truncation [in [§], we obtained M(0"^+)/e^ ~ 2.1 for P G [5,8) at third order 
approximation using the connected graphs as an independent basis]. Our 
value for the lightest glueball can be compared with Samuel's recent result 



1^ from the 2+1 D Hamiltonian QCD in the continuum: M(0~^~^)/e ~ 
1.84 ±0.46. 

One may also understand the relation between the glueball mass and the 
confinement scale from the vacuum wavefunction. The vacuum functional. 



which interpolates the strong and weak coupling regimes, are [pi| , 



\Q) = exp{^ J d^-^x tr[J^,,{VkVk + r')"'/%]}, (15) 

with !F being the field strength tensor in spatial dimensions and V the covari- 
ant derivative. The correlation length ^, with dimension of inverse mass, is 
proportional to e~^, i.e., the confinement scale in the vacuum. It is suggested 
in ref. [|12| that might also be related to the constituent gluon mass and 



the lightest glueball mass. In the strong coupling limit or large Nc (number 
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of colors) limit, it reduces to the strong coupling wavefunction obtained by 



im, 



\n) = exp[-fio J d^-^x trj^\ (16) 

In the intermediate and weak coupling, it becomes [21, ^, ^] 

\n) = exp[-^ I d^'-'x trJ'^ - ^ 1 d^-^x triVJ'fl (17) 

which is just our vacuum wavefunction for the long wavelength configurations 
0]. The correlation length has a relation with the coefficients /io and ^2'- 

i = {^fl\ (18) 

For 2+1 D SU(2), i = O.GS/e^ (see refs. [|l|, |), while for 2+1 D SU(3), 
our result |]^ is ,^ = 0.53/e^. If the glueball mass is proportional to the 
constituent gluon mass, from the difference of the scales between SU(2) and 
SU(3), one may also guess M{0^~^)/e'^ ~ 2, consistent with the result (0) 
from our practical calculation. 



Combining the most recent Monte Carlo data for the string tension 
a in QCDg, which is y/a = (0.554 ± 0.004)e^, we obtain the ratio of the 0++ 
glueball mass over square root of the string tension in the continuum limit 

Mo+^ 



3.88 + 0.11. (19) 



Now we follow the argument of dimensional reduction [|1^, [11], |T2[. In a 
confining theory in D space-time dimensions with 2 < D < 4, the function 
in the exponential of the vacuum functional acts like an action of an effective 
field theory in D — 1 dimensions. In other words, in computing vacuum 
expectation values, a confining theory in D dimensions becomes a localized 
field theory in d = D — 1 dimensions. This can be exactly proven in the 
strong coupling or large Nc limit. Because in this limit, the fixed time vacuum 
expectation value of a operator 0{U) in D dimensions is 

(^IO(^)I^) ^ ndU]0{U)exp[-2f,oJd^-'x trT\ 
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corresponding to the path integral expression for < 0{U) > in D — 1 di- 
mensional lattice field theory. It has been argued [|l^, |2^, |25| that 3+1 D 
theory can still be approximated by its 2+1 D theory for long wavelength 
configurations in comparison to the confinement scale. According to this ar- 
gument, M{J^'-^)/ ^/a for the lightest glueball should be approximately the 
same for 2+1 and 3+1 dimensions. Since for SU(3) the number of color is 
larger, and the measured length ^ in the vacuum functional ([T7|) is smaller 
than that for SU(2), our speculation is that the approximation is better 
for SU(3) gauge theory. In fact, equation (|^) is consistent with IBM data 
M{0++)/^a = 3.95 from Monte Carlo simulation of 3+1 D lattice QCD, 
providing ^/cr = 0.44 Gev is used. From this world average value for the 
string tension and (]TP|), we expect 

M(0++) = 1.71 ± 0.05 GeV, (21) 

in nice agreement with the IBM data M(0++) = 1.740 ± 0.071 @]. This 
favors ^//j(1710) candidate of the 0"*"+ glueball. 

In conclusion, using the eigenvalue equation method developed in P, 0], 
we obtain accurate result for the lightest glueball mass in (2+l)-dimensional 
QCD, with systematic uncertainty under well control. We also use the idea 
dimensional reduction to extrapolate the results to QCD in 3+1 dimensions. 
Perfect agreement with the most accurate Monte Carlo data indicates that 
QCD3 is not just a toy model for QCD4, and their relation should be more 
deeply understood. 
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